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Context. Rotation is known to have a strong impact on the nucleosynthesis of light elements in massive stars, mainly by inducing 
mixing in radiative zones. In particular, rotation boosts the primary nitrogen production, and models of rotating stars are able to 
reproduce the nitrogen observed in low-metallicity halo stars. 

Aims. Here we present the first grid of stellar models for rotating massive stars at low metallicity, where a full s-process network is 
used to study the impact of rotation-induced mixing on the neutron capture nucleosynthesis of heavy elements. 
Methods. We used the Geneva stellar evolution code that includes an enlarged reaction network with nuclear species up to bismuth to 
calculate 25 M models at three different metallicities (Z = 1CT 3 , 10~ 5 , and 1(T 7 ) and with different initial rotation rates. 
Results. First, we confirm that rotation-induced mixing (shear) between the convective H-shell and He-core leads to a large production 
of primary 22 Ne (0.1 to 1% in mass fraction), which is the main neutron source for the s process in massive stars. Therefore rotation 
boosts the s process in massive stars at all metallicities. Second, the neutron-to-seed ratio increases with decreasing Z in models 
including rotation, which leads to the complete consumption of all iron seeds at metallicities below Z = 10~ 3 by the end of core 
He-burning. Thus at low Z, the iron seeds are the main limitation for this boosted s process. Third, as the metallicity decreases, the 
production of elements up to the Ba peak increases at the expense of the elements of the Sr peak. We studied the impact of the 
initial rotation rate and of the highly uncertain [1 0(a, y) rate (which strongly affects the strength of 16 as a neutron poison) on our 
results. This study shows that rotating models can produce significant amounts of elements up to Ba over a wide range of Z, which 
has important consequences for our understanding of the formation of these elements in low-metallicity environments like the halo of 
our galaxy and globular clusters. Fourth, compared to the He-core, the primary 22 Ne production induced by rotation in the He-shell 
is even higher (greater than 1% in mass fraction at all metallicities), which could open the door for an explosive neutron capture 
nucleosynthesis in the He-shell, with a primary neutron source. 
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1. Introduction 

Massive stars in the early universe are different from the ones 
observed today. Their low metal content makes them more com- 
pact and rotate faster than their equivalents found in the Milky 
Way. This view is supported by observations of an increasing 



2007) and by more rapidly rotating, massive stars in the SMC 
compared to the Milky Way (Hunter et al. 2008). Rapidly ro- 
tating stellar models at low Z have been calculated by |Meynet 
|et al.| ( p006l > and |Hirschi ( 2007| l. In these models, primary nitro- 
gen yields are much larger than in non-rotating models. When 
yields from these rotating models are used as input in chemi- 
cal evolution models, a nice fit of the N/O in very metal poor 
halo stars (see e.g. |Spite et al. 2005 \ can be obtained ( |Chia ppini 
et al. 2006 1. This provides a strong support for the occurrence 
of rotation-induced mixing at low Z and for our models. The 
primary nitrogen production in rotating low-Z stellar models is 
accompanied by the primary production of other isotopes like 
13 C, and especially 22 Ne, which is the neutron source for s pro- 
cess in massive stars (e.g. Kappeler et al. 201 1| and references 
therein). 

A first attempt to assess the impact of rotation on the s- 
process nucleosynthesis in low-Z massive rotating stars was 
made by Pignatari et al. ( 2008 1, who investigated the impact of 



large s-process enhancements in one of the oldest globular clus- 
ters in the bulge of our galaxy supports the view that massive 
stars could indeed also be important sources for these elements 
(Chiappini et al. 201 1 1, highlighting the need for comprehensive 
calculations of the s process in low-Z massive rotating stars. In 
this letter, we present the first such low-Z, massive rotating star 
models including a full s-process network. 

We describe our models in Sect. [2] The impact of rotation on 
the s process in massive 25 M Q stars at different metallicities is 
discussed in Sect. [3] We give our conclusions in Sect. |4] 

2. Stellar evolution models 

We calculated our models with the Geneva stellar evolution code 
(GENEC), which is described in detail in (Eggenberger et al. 



2008 ). The main improvement brought to GENEC for this work 
is the integration of a large nuclear reaction network (613 iso- 
topes up to the end of He-burning and 737 from there on). The 
smaller network is almost identical to the s-process network 
used by |The et alT] ( p000] see their Table 1). GENEC with the 



primary Z/ Ne in a parametrised way. The recent observation of 



enhanced nucleosynthesis network was first used to study the 
evolution of the abundance of light elements at the surface of 
main sequence massive stars (Frischknecht et al. |2010 1. Since 
rotation-induced mixing is prime importance in this work, we 
briefly review here the input physics used. We used the horizon- 
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tal diffusion coe fficient of |Zahn (|1992| l and the shear diffusion 
coefficient from |Talon & Zahnf i 1997 ), which is a conservative 
choice since this prescription includes a strong reduction of mix- 
ing across mean molecular weight gradients. 



We used the reaction library (REACLIB) from Rauscher & 
Thielemann (2000), but with major updates. The charged parti- 
cle reaction rates from |Angulo et al.| ( |1999| l were used except 
for the following reactions: /z Ne(a;,n) was taken from Jaeger 
et~aT1p00T] i and the 3a-rate from |Fynbo et al.|p005) . If avail- 
able, the neutron captures were taken from the KADoNiS com- 
pilation ( Dill mann et al.| |2006]> and the temperature-dependent 
beta decays from |Takahashi & Yokoi| ( |1987j i. Several reac- 
tion fits were downloaded from the JINA-REACLIB website 
(groups. nscl.msu.edu/jina/reaclib/db). 

For 17 0(a, y) and 17 0(or, n) reaction rates, we used the rates 
of|Caughlan & Fowler| ( [l9M] l (hereafter CF88) and |Angulo etaL 
( 1999| l, respectively. Their ratio determines the strength of 10 O 
as a neutron poison a nd is very uncertain at the moment. Indeed, 
Descouvemont| ( |1993| predicts that the ll O(a, y) should be a fac- 
tor of 1000 smaller than the CF88 rate. This huge unc ertainty 
strongly affe cts the s process in massive stars at low 2Q(Hirschi 



et al. 



12008b, where ls O is known to be a strong neutron ab- 
sorber/poison (Rayet & Hashimoto] 2000 1. Rece ntly, two inde- 
pendent groups have measured the Q(a, y) rate (|Taggart et al 



201 1 Best et al. 201 l| l, but it is not yet clear whether the new 
rate will be lower than the CF88-rate at the relevant energies 
(priv. comm. A. Laird). To assess the impact of a decreased rate, 
we also calculated 25 M models with the CF88 rate divided by 
a factor 10, which will probably still be inside the uncertainties 
of the new measurements in the relevant energy range. 

We calculated models of 25 M Q stars at initial metallicities, 
Z = 10~ 3 , 10~ 5 , and 10~ 7 (the corresponding [Fe/H]-values are 
given in Table [T}, to investigate the influence of rotation on the 
s process in massive stars. At all three metallicities, we assumed 
an a-enhanced composition^] All other elements were scaled 
from the solar composition published in Asplu ncTet al.| ( [2003] l. 
At each Z, a model without rotation was calculated for reference. 
For rotating models, we chose v m i/u a i t = 0.4 as a standard rota- 
tion rate, where i/; n ; is the initial surface velocity and v crlt the crit- 
ical rotation velocity, for which the centrifugal force balances the 
gravity at the equator. This ratio corresponds to an average equa- 
torial velocity of (v eq )ivis ~ 214 km s _1 on the main sequence at 
Z Q in our models. This is slightly lower than the peak velocity 
of the velocity distribution of B-type stars in our Galaxy found 
by Dufton Tt al.| ( |2006) . A constant i/mi/fcrit leads to a slight de- 
crease in initial angular momentum with decreasing Z and in our 
models at lowest Z to (f e q)MS ~ 383 km s _1 . At low Z, massive 
stars could have a higher initial angular momentum according to 
the latest star formation simulations (Stacy et al. 2011 
also calculated models with i/ini/fcrft = 0.5 for Z 
fini/fcrit = 0.6 for Z = 10~ 7 , corresponding to (v eq ) M s 
and 588 km s _1 , respectively. All models were computed un- 
til at least the start of O-burning. Beyond this point, the carbon 
and helium-burning shells do not change significantly, so the s- 
process yields can be estimated well at this point. The yields can 
still be slightly modified by late interactions between shells and 



1 Several other rate uncertainties e.g. 22 Ne(or, y) and 22 Ne(or, n) also 
affect the results and will be studied in a forthcoming paper. 

2 with the a-elements enhanced with respect to iron, i.e. [X/Fe] = 
A = constant (since [Fe/H] < -1) where A = +0.562, +0.886, +0.500, 
+0.41 1.+0.307, +0.435, +0.300, +0.222, and +0.25 1 for I2 C, 16 0, 20 Ne, 
24 Mg, 28 Si, 32 S, 36 Ar, 40 Ca, and 48 Ti, respectively 



so we 
10~ 5 and 
428 



increased temperature just before the collapse and by explosive 
nucleosynthesis (see e.g. Tur et al. 2009 1. 

We list the properties of the models in Table [TJ in the fol- 
lowing order: model name (row 1), initial metallicity (2), [Fe/H] 
(3), i^ini/fcrit (4), average rotation rate on the main sequence 
(v)ms (5), the mass fraction of 22 Ne burned during core He- 
burning AX( 22 Ne) c (6), the mass fraction of 22 Ne at the start 
of C-burning, X( 22 Ne) (7), the average number of neutron cap- 
tures per iron seed n c (8), and the pre-supernova yields in solar 
masses and production factors, f, for the elements Sr (9/10), and 
Ba (1 1/12). The remnant masses used to calculate the yields (see 
method in Hirschi 2007) range between 2.3 and 2.8 M Q and are 
usually situated close to the bottom of the last carbon shell. 

We can see from Table [l](AX( 22 Ne)) that rotating models at 
all metallicities produce and burn significant amounts of 22 Ne, 
confirming the results of previous studies ( |Meynet et al.| 2006 



Hirschi| [2007 ; Pignatari et al.| [2008 ). The amount of primary 
Z/ Ne in the convective He-core at the end of He-burning when 
s process is active is between 0.1 and 1% in mass fractions. 
Considering a constant value of v m i/u a -i t = 0.4 at all metal- 
licities (models Al, Bl, CI), the primary 22 Ne in the He-core 
decreases slightly with decreasing metallicity. There is, how- 
ever, theoretical and observational support to consider a slight in- 
crease in i/ini/i'crit with decreasing metallicity. Considering mod- 
els Al (Ufoi/i/ait = 0.4 at Z = 10~ 3 ), B3 (w ini /v crit = 0.5 at 
Z = 10~ 5 ), and C3 (Vmi/fcrit = 0.6 at Z = 10~ 7 ), which corre- 
spond to an increase in fmi/fcrit with decreasing metallicity, we 
see that rotating models produce and bum a constant quantity 
of 22 Ne, around 0.5% in mass fraction. We can see from these 
results that rotating models produce significant amounts of 22 Ne 
in a primary way over the entire range of metallicity computed, 
whereas in non-rotating models (AO, B0), the amount of 22 Ne 
available is secondary. 

In the He-shell, the 22 Ne mass fraction in rotating models is, 
by a factor of 19 (Al), 1100 (Bl) and 13400 (CI) higher than 
the corresponding non-rotating models. While the 22 Ne in the 
He-shell of non-rotating model is purely secondary, in rotating 
models it is primary and amounts to 1 and 2% in mass fraction at 
the pre-SN stage and independent of metallicity. This is very in- 
teresting for explosive nucleosynthesis in the He-shell. This site 



was investigated for the r-process by Truran et al. ( 1978| ) and 
Thielemann et al. ( 1979| ), but later on considered only as a pos- 



sible site for a weak r-p rocess ( |Woosley et al.||2002] ) or n-process 
(Rauscher et al. 2002}. Given the large amounts of primary 22 Ne 



produced in the rotating models at all Z, it is worthwhile recon- 
sidering explosive nucleosynthesis in the He-shell. 

3. Non-standard s process 

The standard s process in massive stars is a secondary process 
(see e.g. Raiteri et al. 1992), because the main neutron source 
( 22 Ne coming from the initial C, N, and O) and the seeds (mainly 
iron) have a secondary origin, while neutron poisons are a mix- 
ture of primary (mainly 16 for He-buming and 2l) Ne, 24 Mg for 
C-burning) and secondary (mainly 25 Mg for He-burning) ele- 
ments. At low Z, the contribution from the carbon shell becomes 
smaller (given the mass fraction of 22 Ne left at the start of C- 
burning, X( 22 Ne), see line 7 in Table]]}. In non-rotating models 
(AO, B0), the standard s process is even less efficient than a sec- 
ondary process. 

Rotation significantly changes the structure and pre-SN evo- 
lution of massive stars ( |Hirschi et al.[ |2004| ) and thus also the 
s-process production. Rotating stars have similar central proper- 
ties to more massive non-rotating stars. In particular they have 
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Table 1. Parameters of the 25 M Q models calculated. 



model" 


AO 


Al 


BO 


Bl 


B2* 


B3 


B4* 


CI 


C3 


C4* 


Z 


10- 3 


io- 3 


10" 5 


io- 5 


io- 5 


io- 5 


IO" 5 


io- 7 


io- 7 


io- 7 


[Fe/H] 


-1.8 


-1.8 


-3.8 


-3.8 


-3.8 


-3.8 


-3.8 


-5.8 


-5.8 


-5.8 


^ini / ^crit 





0.4 





0.4 


0.4 


0.5 


0.5 


0.4 


0.6 


0.6 


{v) h d 





285 





333 


333 


428 


428 


383 


588 


588 


AX( Ne) 


7.52(-4) 


4.08(-3) 


7.21(-6) 


1.23(-3) 


1.27(-3) 


3.83(-3) 


3.75(-3) 


1.05(-4) 


4.57(-3) 


4.44(-3) 


X( 22 Ne) 


4.16(-4) 


6.22(-4) 


4.14(-6) 


1.69(-4) 


1.82(-4) 


4.94(-4) 


4.85(-4) 


3.81(-5) 


2.68(-4) 


3.11(-4) 


< 


1.64 


12.7 


0.08 


5.77 


23.1 


16.5 


31.8 


14.0 


33.5 


48.5 


pre-SN yields [M Q ] 


Sr 


1.66(-8) 


1.13(-5) 


2.60(-ll) 


8.70(-9) 


3.82(-7) 


2.05(-7) 


4.36(-7) 


1.08(-9) 


3.66(-9) 


3.02(-9) 


f(Sr) 


1.75 


501.4 


1.12 


28.86 


1724 


944.7 


1968 


489.1 


1650 


1361 


Ba 


2.13(-9) 


2.29(-8) 


3.56(-12) 


7.41(-11) 


6.49(-9) 


4.89(-10) 


4.73(-8) 


3.80(-10) 


9.86(-10) 


2.95(-9) 


f(Ba) 


1.34 


3.87 


1.06 


2.10 


106.1 


9.07 


767.0 


615.6 


1598 


4771 



Notes. <fl) Model names (1 st line) include a letter and a number. The letter refers to the initial metallicity, Z (2 line) and the number to the initial 
rotation/170(or, y) rate used. is for non-rotating models, 1 for models with constant v-^jv^ = 0.4, 3 for increasing v^dv<. Ax : 0.5-0.6 with 
decreasing Z and 2 and 4 for lower 170(a,y) rates (CF88/10). ih) Model with the \10{a,y) reaction rate of CF88 divided by 10. w in km s" 1 . (d) 
Since 22 Ne is produced and destroyed at the same time in rotating stars, we derived this value from the sum of the 2, Mg and 26 Mg produced during 
central He-buming. (c> Number of neutron captures per initial seed. 



more massive helium-burning cores and higher central tempera- 
ture. The latter means that the He-core contribution increases at 
the expense of the C-shell contribution. This and the primary ori- 
gin of neutron poisons during carbon burning lead to a very small 
C-shell contribution to the pre-SN yields at low Z, i.e. lower than 
10% for all the rotating models presented in Table [T] 

All the effects of rotation boost the s-process production at 
all metallicities compared to non-rotating models as can be seen 
in Table[l] At [Fe/H] = -1.8 (models AO and Al), the rotating 
model burns about five times more 22 Ne than the non-rotating 
model during core He-burning. Although a factor of 5 is not very 
large, 96% of the initial iron is transformed into s-process ele- 
ments in the rotating model, whereas only 53% of the initial iron 
is transformed in the non-rotating model. In the rotating model, 
the neutron source is sufficient to deplete all the seeds, and the 
production is limited by the seeds (not the neutron source any 
more). This explains why the s-process yields are much larger 
in the rotating model, by up to a factor of 1000. Nevertheless, 
at [Fe/H] = -1.8, rotating models produce elements in large 
quantities only up to the Sr peak. We thus expect that the non- 
standard s process is qualitatively similar to the normal s process 
(production of elements with A = 60-90) for [Fe/H] > -2 but 
with higher production factors. At [Fe/H] = -3.8 (see Fig. [T] 
showing production factors at the pre-SN stage), elements up to 
Sr are still strongly produced in the rotating model (Bl, Fig. [T]l, 
whereas the overproduction is very small in the non-rotating 
model. Depending on the initial rotation rate and the n O(a, y) 
rate used, elements up to Ba start to be efficiently produced (see 
discussion below). At [Fe/H] = -5.8, both the rotating and non- 
rotating models produce and burn enough 22 Ne to deplete all the 
iron seeds. At this very low Z, the main limitation is the seeds 
and, given the low initial iron content, the s-process production 
remains very modest. 

It is interesting to look at the metallicity dependence of the 
non-standard s-process production in rotating models. As said, 
the production is limited mainly by the iron seeds. Even at the 
lowest metallicities in a very rapidly rotating model (C3 and 
C4, Vi„i/iA;rit= 0.6 instead of the standard 0.4), and thus with 
a larger primary neutron source, there is no additional produc- 
tion of s-process elements starting from light element seeds like 
22 Ne. Instead, what happens is that not only iron is depleted, but 
also elements up to Sr are partially destroyed (after being pro- 



duced) and heavier elements like Ba are produced ([Sr/Ba] < 0). 
Indeed, going from [Fe/H] = -3.8 (Bl) to [Fe/H] = -5.8 (CI), 
the Sr yield decreases by a factor of » 9, while the Ba yield 
increases by a factor of 5. We therefore have a different metal- 
licity dependence for the production of elements belonging to 
the different peaks: there is a roughly secondary production of 
elements up to Sr but the Z dependence for heavier elements 
like Ba is milder. The secondary-like behaviour of Sr/Y/Zr in 
the metallicity range covered by our models, makes the non- 
standard s process in massive rotating stars an unlikely solution 
for the LEPP (light element primary process) problem at low Z 
( |Travaglio et aLl|2004| . 

Apart from the seeds, the s process strongly depends on the 
neutron source and neutron poisons. Both are still quite uncer- 
tain at present. The neutron source, 22 Ne, depends on rotation- 
induced mixing, hence on the initial velocity of stars at very low 
Z. To assess the uncertainty linked to the primary production of 
22 Ne, we calculated additional models at [Fe/H] = -3.8 and 
-5.8 with a higher initial velocity (B3/B4 v^/v^ = 0.5 and 
C3/C4, Uini/wcrit = 0.6). Compared to the CI case (fini/i'crit= 
0.4), the increase in the production of 22 Ne is about a factor of 
4. This leads to a higher neutron capture per seed (n c ) and thus 
to a production of elements like Ba at the expense of elements 
like Ge, but the total production (sum of all isotopes heavier than 
iron) is still limited by the iron seeds, as said above. A major un- 
certainty concerning neutron poisons is the importance of 16 as 
a neutron poison. As explained in the previous section, at low 
Z, ls O is a strong neutron absorber during core He-burning. The 
neutrons captured by 16 0(n, y) 17 may either be recycled via 
17 O(a,n) 20 Ne or lost via 17 0(a,y) 21 Ne. To investigate the im- 
pact of the I7 0(a,y) 21 Ne uncertainty, we calculated models with 
the CF88 rate divided by a factor of 10 (B2, B4, and C4). The 
impact of a change of even a factor of 10 in this rate is strong 
(see Table [T] and Fig. [TJ. Given the differences between models 
with the CF88 and CF88/10 rate, the experimental determination 
of the n O(a, y)-rate and the 17 0(a, n) is crucial to a more accu- 
rate prediction for the s process in massive rotating stars at low 
metallicity. 

Characteristic quantities for the s process are [Sr/Ba] and 
[Sr/Pb]. For model Al at Z = 10~ 3 , [Sr/Ba] = +2.33, which 
is typical of the standard s process in massive stars. On the other 
hand, the models at lower metallicities show reduced values be- 
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Fig. 1. Production factors (ejected mass divided by the initial mass of the element) for the 25 M Q models with Z = 10~ 5 . The model 
without rotation (triangles) does not produce s-process efficiently, whereas the rotating models (circles, Bl, and diamonds, B3) do. 
The additional, rotating models with reduced ll O(a,y) rates (stars, B4, CF88/10) highlight the uncertainty for the neutron poison 
"'(). 



tween +2.05 and +0.42 at Z = 10~ 5 , and between +0.02 and 
-0.54 at Z = 10~ 7 (although the Z = 10~ 7 models will not con- 
tribute much to the galactic chemical enrichment). We therefore 
generally obtain [Sr/Ba] > from the low-Z, massive rotating 
star, s process. Pb is only produced in very extreme cases and 
with [Pb/Sr] < -1. Thus massive stars cannot be responsible for 
CEMP-s stars. The models of Pignatar l et al.| p008) efficiently 
produce elements up to Pb. This difference between the two stud- 
ies is explai ned by the fact that [Pignatari et al.| ( |2008[ ) used the 
ll O(a, y) of Descouvemont ( 1993 1, which is now disfavoured 
by recent experiments. Finally, a similar spread for [Sr/Ba] is 
obtained as the one used in Chiap pini et al.| ( [201 \) to explain the 
abundance of one of the oldest globular clusters in the galactic 
bulge, thus supporting the conclusions of that study. 

4. Conclusions 

We have calculated the first complete stellar models including 
a full s-process network to study the effects of rotation on the 
s process in massive stars in low-metallicity environments. Our 
models confirm that rotation-induced mixing in radiative zones 
leads to a primary production of 22 Ne in the He core over the 
wide Z range studied ([Fe/H] = -5.8, -3.8, and -1.8). From 
0.1 to 1% of 22 Ne in mass fraction is produced in central He- 
buming. The large primary production of 22 Ne also occurs in the 
He-shell in all the rotating models calculated, i.e. 1-2% in mass 
fraction at all Z. This can have a strong impact on the super- 
nova nucleosynthesis in He-burning shell, since it could provide 
a primary neutron source for neutron capture nucleosynthesis at 
all values of Z. We investigated the impact of the ll O(a, y)-rate 
on the strength of 16 as a neutron poison and find the impact 
of current uncertainties still significant, thus deserving further 
experimental studies. 

The s process in massive stars with rotation contributes sig- 
nificantly to the production of elements up to Sr above [Fe/H] w 
-2, the main limiting factor being the iron seeds at low Z. The 
production of elements up to Ba is sensitive to the initial ro- 
tation rate, and we generally obtain [Sr/Ba] > from low-Z 
massive rotating star s process. Pb is only produced in very ex- 
treme cases, but with [Pb/Sr] < -1 massive stars cannot be re- 
sponsible for CEMP-s stars. A similar spread for [Sr/Ba] to the 
one used in Chiapp ini et al.| ( |201 l| l was obtained to explain the 



unique abundance of several stars in one of the oldest globular 
clusters in the galactic bulge, thus supporting the conclusions of 
that study. Further stellar and galactic chemical evolution mod- 
els will assess the full impact of this boosted s process. 
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